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Abstract 



The theory of preasymptotic effects in inclusive decays of heavy flavors is briefly 
reviewed. 



1 Introduction 



Heavy flavor hadrons Hq contain a heavy quark Q plus a cloud built from light 
quarks (antiquarks) and gluons. The heavy quark Q experiences a weak transition. 
The nature of this transition is of no concern to me here. It can be a radiative 
transition, like b — > sj, semileptonic decay like b — > clvi or a non-leptonic decay to 
lighter quarks, e.g. b — > cud. It is assumed that at short distances the amplitude 
is known from the electroweak theory. The task of the QCD-based theory is to 
calculate preasymptotic effects in the decay rate of the hadron Hq and other decay 
characteristics: the energy spectra, the average invariant mass of the hadronic state 
produced, etc. These effects are due to interactions with the soft degrees of freedom 
in the light cloud. 

The foundation of the theory was laid in the eighties [|l[] when it was realized 
that the operator product expansion could be used in application to the so called 
transition operator of the type 

f(Q^ f - Q) = z J d'xiCwix), £^(0)} T , (1) 

where Cw is the short-distance weak Lagrangian governing the transition Q — > / 
under consideration. The momentum operator of the heavy quark Q is written 
as a sum of two terms, P M = mqv^ + tt^ where t> M is the four- velocity of the heavy 
hadron Hq, 7r M is the residual momentum operator, responsible for the interaction 
with the "background" gluon field in the light cloud. The large mechanical part 
vrtQV^ in the heavy quark momentum guarantees that the transition operator ([!]) 
can be found as a sum of local operators (with some reservations to be discussed 
below). These operators are ordered according to their dimensions. At the level of 
l/mn we have only two operators; extra four-fermion operators are added at the 
level of I/tuq. A few of these were calculated 8 years ago 

At the next stage each term in the expansion must be averaged over the hadronic 
state Hq. At this stage the bound state dynamics is accounted for. 

After the initial excitement the OPE-based theory of the inclusive heavy flavor 
decays was in a rather dormant state until recently. The revival it experiences now 
is due to a combination of several factors. First, a very concise and convenient 
language was created, the heavy quark effective theory (HQET) j3|. Calculations of 
the non-perturbative effects were translated in this language and developed in Refs. 
||, ||, ||. Second, all relevant operators appearing at the level up to 1/mg were 
catalogued and our understanding of their matrix elements (I mean the numerical 
values) was significantly advanced. Finally, the issue of convergence of the non- 
perturbative series was clarified. In the beauty family one expects that the first 
two or three terms in the expansion ensure reasonable accuracy of the predicted 
lifetimes. As for the charmed quark we will see that it is, perhaps, too light for 
duality to set in. Since OPE is used in the Minkowski domain the validity of duality 
is crucial for the whole approach. Theoretically the onset of duality is correlated 
with the behavior of high-order terms in the non-perturbative series. 



1 



2 Master equation 



The OPE-based approach is applicable in a very wide range of problems. Here we 
will concentrate on the total inclusive widths. Generically the tuq 1 expansion for 
the width has the form (for defmiteness I will speak about the beauty family) 



H 2 m 5 

v JJ 1927T 3 ' 1 



X 



f AH b \bb\H b ) C5 (/) (H b \b(i/2)cGb\H b ) 
1 C3UJ 2M Hb + ml 2M Hb 

^ cfjf) (H.KbT^jqT^lH,} + 0( ^_ 4) 1 _ (2) 

The coefficient functions Cj(/) depend on the particular inclusive transition con- 
sidered and are calculable. They are determined by short- distance QCD provided 
that the energy release is large enough. On the other hand, the matrix elements on 
the rhs describe the response of the soft degrees of freedom on the instantaneous 
perturbation, the 6-quark decay. These quantities are essentially non-perturbative. 
But they are universal, and, as seen from Eq. (fj), there are only a few of them. 

The matrix element of the chromomagnetic operator oG is expressible through 
spin splittings, say Mb* — Mb- The four-fermion operators of dimension 6 can 
be evaluated, in the case of mesons, within factorization. For baryons a reliable 
calculation of the corresponding matrix elements is a problem essentially unsolved 
so far. As for the scalar density, 66, it is this term that exactly reproduces the parton 
model (asymptotic) result, plus preasymptotic corrections. This operator also can 
be written as an expansion, 

66 = 6 7o 6 \rb(n 2 - (i/2)aG)b+ 

-L r g 2 b l0 T a bJ2™ Ta <i + 0(l/mt) (3) 

4m Q g 

where the sum runs over the light quarks. A new operator appearing here is 67? 2 6, 
the square of the spatial momentum of the b-quark. The matrix element of this 
operator in the B meson can be limited from below, for a detailed discussion see 
Ref. fTj. The average spatial momentum turns out to be surprisingly large, larger 
than 0.6 GeV! The QCD sum rule calculations [H] yield even a larger value, ~ 0.7 
GeV. The expectation value of 7? 2 in baryons is expected to be close to that in 
mesons. 

Time/space limitations do not allow me to go into further details. Let me point 
out only the most remarkable features of the overall picture. 

(i) The total rates do not contain non-perturbative corrections of order l/m b , the 
so called CGG/BUV theorem 0, ||. The corrections start at the level This 
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sets the scale of preasymptotic effects in the beauty family at the level of several per 
cent since (B\biaGb\B) /2m\ ~ 0.03. In particular, deviations of the lifetime ratios 
from unity are expected to be of this order of magnitude. At the level l/mf all B 
mesons have the same lifetimes (disregarding some small SU (3)a breaking effects). 

(ii) The difference in the lifetimes of baryons and mesons is due to the fact that 
the expectation values of the operators in Eq. (|2]) are different for mesons and 
baryons. This difference arises at the level 

(iii) Four-fermion operators of dimension 6 produce effects formally scaling like 
1/mj- , although numerically they seem to be enhanced since the corresponding coef- 
ficients have one loop less and, additionally, a key constant fs turns out to be rather 
large. This enhancement may lead to the fact that dimension 5 and 6 operators are 
competitive in the beauty family. The four-fermion operators shift the lifetimes of 
mesons versus baryons and split the meson lifetimes from each other. 

(iv) Situation with B s is exceptional. The lifetime difference between -B SjS hort and 
-Bs,iong is due to a mechanism not exhibited in Eq. (Q), namely B — B oscillations. 
The corresponding estimates were done in Ref. §|. 



3 Phenomenological implications 

Assembling all theoretical elements discussed above (and those which are discussed 
in the original literature) we arrive at the following pattern. The lifetime of a charged 
B meson is predicted to exceed that of a neutral B meson, 

-j^l ~ 1 « 0.05(/ B /200MeV) 2 ~ 0.05. (4) 

At this level it is expected that r(fi^) w f(B s ) where f denotes the average lifetime 
of the two mass eigenstates in the B°-B° system. It is curious that B s oscillations 
will seemingly produce the largest lifetime difference, 

^§4- - 0.18(/ B /200MeV) 2 ~ 0.18 . (5) 

The baryon matrix elements are always most difficult for consistent analysis; there- 
fore, the baryon-to-meson lifetime ratios should be taken with caution. Still, plau- 
sible estimates indicate that one can expect T(Ab)/r(Bd) ~ 0.9. 



4 A grain of salt: Bt s \(B) 

The theory of preasymptotic effects which I have just sketched, being applied to the 
problem of the semileptonic branching ratio in the B mesons, leads to a paradox. 
In this case the heavy quark expansion can be readily carried out up to terms 
of order 1/mf. One obtains a formula very similar to Eq. (|j), with the same 
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structure and the same operators [I0|. The leading non-perturbative correction 
0{m^ 2 ) tends to diminish the branching ratio while the term 0(m^ 3 ) tends to 
increase it. Both effects, however, are far too small to produce a noticeable impact 
on the branching ratio. At best they shift the prediction for the branching ratio 
by 0.5% or less. Thus we are forced to conclude that the prediction for Br s \(B) is 
controlled by perturbative QCD. People believe that perturbative QCD typically 
yields Br sl (£>) « 13%; twisting arms allows one to go down to 12.5% [p| . At 
the same time experimentalists, both CLEO and ARGUS, seem to be firm in their 
conclusion that Br s \(B) < 11%. A natural question is what went wrong? 

I leave aside the possibility that the experimental numbers are wrong. There is 
no visible loophole in the OPE-based theory of preasymptotic effects either. Then 
the remaining logical options are as follows: (i) something is missing in the pertur- 
bative analysis; (ii) new physics shows up in the B meson decays. Both options 
must be investigated. In a recent paper a new contribution in the perturbative 



calculation is identified, not included in the analysis of Ref. [fLl], which seemingly 
works in the right direction - diminishes Br s i(i?)p er t by ~ 0.5%. It remains to be 
seen whether the perturbative number can reach the 11% mark under realistic choice 
of relevant theoretical parameters (the quark masses, a s , etc.). (Let me make a side 
remark: I do not believe that a s (Mz) can be as large as 0.126, as is allegedly implied 
by the so called global fits at the Z peak at present. A wealth of low-energy data 
point to a significantly lower value of a s , something like 0.114 or even lower. In terms 
of Aqcd the difference is drastic. I urge to take this discrepancy very seriously.) 



5 The family of charm 

It might seem to be a trivial exercise to substitute by m c in the master equation. 
Yes, technically this is easy, and formally all 1/m^ and \jvn? c corrections have been 
written down and classified. They are much larger, of course, than in the beauty 
family; typically of order of 0.5. I refer those interested to a very detailed recent 



update [13]. Qualitatively the pattern of the lifetimes in the charm family emerging 
in the heavy quark expansion agrees with experiment. Namely, those particles that 
live longer are predicted to live longer, etc. However, quantitatively the 0{m~ 2 ) and 
0(m~ 3 ) preasymptotic terms are smaller than what one needs in order to reproduce, 
say, r(D + )/r(S°) cxp ~ 12. Since arithmetically the calculation is certainly correct 
one may start suspecting that something went wrong in the basics. 

The operator product expansion, the foundation of the whole approach, is a 
well defined procedure in the Euclidean domain. A specific feature of the transition 
operator (^ is its essentially Minkowski character. Therefore, in justifying the short- 
distance calculation of the coefficient functions one must always keep in mind a kind 
of analytic continuation, through a dispersion relation. Thus, strictly speaking, 
theoretical predictions for Cj(/) in Eq. (|2|) refer to quantities integrated over energy 
in some energy range. 
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If the energy release is large enough so that duality is valid and the integrand is 
smooth, this smearing is unimportant; one can predict the coefficient functions for 
the given energy release, locally. It is always tacitly assumed that this is the case. 
The onset of duality is governed by exponential terms, not visible to any finite order 
in 1/m.Q expansion. Due to this reason the onset must be abrupt. 

We are inclined to think that the c quark is not heavy enough to warrant dual- 
ity. The strongest argument comes from consideration of T S \(D). Indeed, with the 
reasonable value of m c (m c (m c ) ~ 1.3 GeV) the parton-model prediction is close to 
the experimental number. The first perturbative correction is negative |l lj and the 
second seems to be negative as well [14j] . The non-perturbative corrections follow the 
same pattern. The leading term is known from Ref. || while the l/m? c cor- 

rection has been estimated recently [T5|]; both are negative. The combined effect of 
the leading corrections amounts to reducing T S \(D) by 50%, and the next-to-leading 
terms only worsen the situation! 

This work was supported in part by DOE under the grant number DE-FG02- 
94ER40823. 
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